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r^ I Abstract: The unfortunate case where the Large Hadron Cohider (LHC) fails to discover 

physics Beyond the Standard Model (BSM) is sometimes referred to as the "Nightmare 
scenario" of particle physics. We study the consequences of this hypothetical scenario 
for Dark Matter (DM), in the framework of the constrained Minimal Supersymmetric 

S^ . Standard Model (cMSSM). We evaluate the surviving regions of the cMSSM parameter 

space after null searches at the LHC, using several different LHC configurations, and study 
the consequences for DM searches with ton-scale direct detectors and the IceCube neutrino 
telescope. We demonstrate that ton-scale direct detection experiments will be able to 
conclusively probe the cMSSM parameter space that would survive null searches at the 
LHC with 100 fb~^ of integrated luminosity at 14TeV. We also demonstrate that IceCube 
(80 strings plus DeepCore) will be able to probe as much as ^ 17% of the currently favoured 
parameter space after 5 years of observation. 



Contents 

|l]. Introduction |l| 

H). Theoretical Setup | 

^ ThecMSSM | 

^^ Statistical framework H 

PI Priors | 

2^ Cosmological and collider constraints ^ 

LHC constraints IBI 



|3|. Results 

^ The cMSSM after null searches at the LHC @ 

Implications for Direct detection S 



Implications for indirect detection with the IceCube neutrino telescope 12 



3.4 The role of the Higgs 17 



|4|. Discussion and Conclusions |T8 



1. Introduction 

One of the main goals of the Large Hadron Collider (LHC) at CERN is to test the existence 
of physics Beyond the Standard Model (BSM), especially in view of the possible connection 
with the problem of Dark Matter (DM) in astrophysics and cosmology [|l| . Recently, the 
ATLAS and CMS experimental collaborations have published the results of the first LHC 
searches for Supersymmetry (SUSY) ||2|, ^, |^, ||, ^, [TJ^ based on ^/s = 7TeV collisions 
recorded during 2010. The absence of any excess above Standard Model (SM) predictions 
allows one to already set interesting constraints on BSM physics, and several groups of 
authors have already studied the impact of these results on Supersymmetric scenarios (see, 

e.g., gig, |n|, Hill, 0). 



The prospects for discovering BSM physics at the LHC and the consequences for Dark 
Matter searches have been thoroughly discussed in the literature. In particular, we have 
recently discussed the case of positive detection both at the LHC and with Dark Matter 



experiments [|^]. Here, we focus instead on the unfortunate case where the LHC fails 
to discover BSM physics, a scenario sometimes referred to as the "Nightmare scenario" 
of particle physics (see, e.g., |[l^ and references therein). We focus on the constrained 
Minimal Supersymmetric Standard Model (cMSSM) 0, |l|, ||, |0|, ^, for which detailed 



studies exist regarding the reach of the LHC for various configurations of beam energy and 
total integrated luminosity. 

The main focus of this paper is to analyse the consequences of the nightmare scenario 
for Dark Matter searches. These can be actually divided into two broad classes, viz. direct 
and indirect searches (for recent reviews see, e.g., [||, 22, |23|, 24, ^). Direct searches are 



based on the search for rare nuclear recoils due to the scattering of Dark Matter particles off 
nuclei in underground experiments. There are many undergoing and upcoming experiments 
(see, e.g., [^ for a discussion of upcoming experimental capabilities), and despite some 
intriguing signals that have been discovered by the DAM A/LIBRA [^ and, more recently, 
the CoGeNT [|^, ^, Q collaborations, it appears difficult to reconcile a standard Dark 
Matter interpretation of these results with other experimental findings, in particular with 
the nuh searches from XENON- 100 Q or CDMS H |3|] (see, e.g., the discussions in 

Indirect searches are based instead on the astrophysical searches for secondary par- 
ticles produced in the annihilation or decay of Dark Matter particles. Despite the huge 
interest attracted by astrophysical data (see, in particular, the plethora of literature de- 
voted to explaining the positron excess recently measured by PAMELA ||3^), a convincing 
identification of Dark Matter based on these data seems problematic, due to the lack of 
strong constraints on the properties of Dark Matter particles as well as to the poor control 
of astrophysical backgrounds and associated systematics (see, e.g., the discussion in [|lq] ). 
Among the possible smoking-gun signals that can be obtained with indirect searches, we 
will focus on one of the most intriguing, and less affected by astrophysical uncertainties: 
the detection of high-energy neutrinos from the centre of the Sun [^ using the IceCube 
neutrino telescope, which is located at the geographical South Pole (see, e.g., ||39[|). 

We show here that upcoming ton-scale direct detection (DD) experiments, and to a 
lesser extent IceCube, can actually probe a large portion of cMSSM parameter space that 
will remain unconstrained in the nightmare scenario, and therefore represent a unique 
opportunity to discover new physics by the end of the decade^ in case of null searches at the 
LHC. 

The paper is organized as follows: in Sec.|| we introduce our theoretical framework, i.e., 
the cMSSM, and we present the statistical tools adopted to identify the region of cMSSM 
parameter space corresponding to the LHC nightmare scenario; in Sec. || we present our 
results and discuss the consequences for Dark Matter searches; and, finally, in Sec. ^ we 
provide a brief summary. 



2. Theoretical Setup 

2.1 The cMSSM 

Supersymmetry (SUSY) is one the most compelling extensions of the Standard Model of 
Particle Physics, and the Supersymmetric neutralino is among the best-motivated Dark 
Matter (DM) candidates (for a review of SUSY DM candidates see, e.g., [|o|, H, gsf.) 

Neutralinos are the mass eigenstate of a mixture of bino, wino (superpartners of the B 
and W^ gauge bosons respectively) and two higgsinos (superpartners of the H^, H2 Higgs 
bosons). In many realizations of the Minimal Supersymmetric Standard Model (MSSM), 
the lightest among the four neutralinos is the lightest SUSY particle, which is made stable 
by virtue of the conservation of R-parity. It is commonly referred to as "the" neutralino, 
Xi, and constitutes a popular DM candidate, adopted by much of the relevant literature. 

Since superpartners of Standard Model particles have not been observed, SUSY must 
be broken at some energy scale. Over the past few decades, several SUSY-breaking mech- 
anisms have been proposed in the literature (for a recent review see, e.g., |Q). Perhaps 
the most appealing of these regards gravity as the means of communication between a 
hidden sector, where SUSY breaking occurs, and the visible sector @. In the so-called 
constrained Minimal Supersymmetric Standard Model (cMSSM) [|^, |l8|, H, |o|, |2l| , which 



is the model focussed upon in this paper, the mechanism via which SUSY is broken nor 
how this breaking is transmitted to the visible sector need be specified, so long as the 
latter occurs at or near the GUT scale. In the cMSSM the soft parameters are assumed 
universal at a high scale (Mx), hence the parameter space of the model is defined in terms 
of four free parameters and one sign: the common scalar mass (mo), the common gaugino 
mass (mi/2) and the tri-linear mass (Aq) parameters (all specified at the GUT scale), plus 
the ratio of Higgs vacuum expectation values tan^S and sign(/i), where fi is the higgsino 
mass parameter whose value is determined from the conditions of radiative electroweak 
symmetry breaking. 

2.2 Statistical framework 

We denote the parameter set of the cMSSM (i.e., mo, mi/2, Aq and tan/3) by 0, where 
we fix sgn(^)=+l, motivated by consistency arguments involving measurements of the 
anomalous muon magnetic moment, which can only be explained with a positive fi [^ (see 
also, e.g., Q), while '0 denotes the relevant SM quantities that enter into the calculation 
of observable quantities (the so-called "nuisance parameters" ) , namely 

^P = {Mt, mhimh)'^, as{Mz)'^, a,^{Mz)'^} , (2.1) 

where Mf is the pole top quark mass, mi){mi))^^ is the bottom quark mass at 7715, while 
aem(Mz)^^ and as{Mz)^^ are the electromagnetic and the strong coupling constants at 
the Z pole mass M^, the last three being evaluated in the MS scheme. We denote the full 
8-dimensional parameter set by 

e = (e,V). (2.2) 
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The cornerstone of Bayesian inference is Bayes' Theorem, which reads 



.im-^-^^. 



The quantity p(Q\d) on the l.h.s. of Eq. (|2.3D is known as the posterior^ while on the r.h.s., 
the quantity p{d\Q) is the likelihood (when taken as a function O for fixed data, d). The 
quantity p{Q) is the prior which encodes our state of knowledge about the values of the 
parameters O before scrutinising the data. Our state of knowledge is then updated to 
the posterior via the likelihood. Finally, the quantity in the denominator of Eq. ( |2.2D is 
known as the evidence or model likelihood. If one is interested in constraining the model's 
parameters, the evidence is merely a normalization constant, independent of O, and can 
therefore be dropped by converting Eq. ( ^ ) into a proportionality relation (see, e.g., [44| 
for further details). 



In order to explore the posterior of Eq. ( |2.3D in an efficient way, we adopt the MultiNest 
algorithm, as implemented in the public package SuperBayeS-vl .5. MultiNest pro- 
vides an extremely efficient sampler even for likelihood functions defined over a parameter 
space of large dimensionality with a very complex structure. This aspect is very important 
for exploring the cMSSM, as previous MCMC scans have revealed that the 8-dimensional 
likelihood surface is very fragmented and that it presents many finely-tuned regions that 
are difficult to explore with conventional MCMC scans (and almost impossible to find with 
conventional grid scans) (see also, e.g., ||4^, E^ for further discussions of these aspects). 



2.3 Priors 

In order to perform our scan over the cMSSM and SM nuisance parameters, we need to 
specify our prior in Eq. (^). The role of the prior is to define a statistical measure on the 
parameter space. In principle, when the likelihood is strongly constraining (i.e., for accurate 
data) the posterior is dominated by the likelihood and the choice of prior is irrelevant, 
as the information in the likelihood completely overrides the information in the prior. 
However, it has been shown that this is presently not the case for the cMSSM, i.e., different, 
plausible, choices of priors lead to different posteriors and hence different inferences on 



cMSSM parameter space [Q, although inclusion of recent direct detection constraints from 
XENON-100 does mitigate the problem [|l4|. It is expected that a complete resolution of 
this issue will come from future, more detailed, data and, in particular, measurements of 
the SUSY mass spectrum by the LHC, which will conclusively resolve ambiguities brought 
about by prior dependences in case SUSY is discovered [^]. 

In this paper, we adopt two sets of widely used priors, namely the so-called "flat priors" 
(uniform on the scalar and gaugino masses) and the "log-priors" (uniform on the log of 
the masses). Both sets of priors are uniform in Aq and tan/3. By considering two choices 
of priors, we can assess the residual prior dependency of our conclusions. The prior ranges 
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Table 1: Experimental data used for the computation of the hkehhood function. For each row, the 
central value is given, together with the experimental and theoretical uncertainty. 



SM (nuisance) 
parameter 


Mean value Uncertainty 
11 a (exp.) 


Ref. 


Mt 

l/aem(M2)^^ 


173.1 GeV 1.3 GeV 
4.20 GeV 0.07 GeV 
0.1176 2 X 10-3 
127.955 0.03 


n 
m 
m 



Table 2: Experimental raean ji and standard deviation <j adopted for the likelihood function for SM 
(nuisance) parameters, assumed to be described by a Gaussian probability distribution function. 



we use are given by 



50 GeV < mo,mi/2 < 4000 GeV 

|Ao|<7TeV 
2 < tan/3 < 62. 
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2.4 Cosmological and collider constraints 

In order to derive current constraints on cMSSM parameters, we perform a global scan 
for both approaches using flat priors on the model parameters, including in the likelihood 
current cosmological and collider experimental constraints, as displayed in Table ||, and 
with constraints on the SM nuisance parameters as displayed in Table|2[ In this paper, we 
utilise the same likelihood function as in [|l4[, invoking the experimental constraints listed 
in Table|^. We refer the reader to |6l| for a detailed description of the likelihood function. 
In particular, we include the measurement of the anomalous magnetic moment of the muon 
based on e^e~ data, which gives a 3.7 a discrepancy with the SM predicted value |5l|] . 

We compute S^^Qj^ at full one-loop level adding the logarithmic piece of the quantum 
electro-dynamics two-loop calculation [^. The BR{B -^ Xgj) branching ratio (which has 
been shown to provide an important constraint, see, e.g., the recent study fO^), has been 
computed with the numerical code SusyBSG |Q using the full NLO QCD contributions, 
including the two-loop calculation of the gluino contributions presented in [^] and the 
results of ||6^ for the remaining non-QCD tan /3-enhanced contributions. For the determi- 
nation of AMbs we use expressions from [|^] which include dominant large tan /3-enhanced 
beyond-LO SUSY contributions from Higgs penguin diagrams. The other B(D)-physics ob- 
servables summarized in Table|| have been computed with the code Superlso (for details 



on the computation of the observables see, e.g., [^ and references therein). Both codes 
have been integrated into SuperBayes. We discard points that do not fulflll the conditions 
of radiative electroweak symmetry breaking and/or give non-physical (i.e., tachyonic) so- 
lutions. 

We also include the constraints on the cold Dark Matter (CDM) relic abundance 



determined from the 7-year WMAP data |57] to constrain the relic abundance Q^h of 



the lightest neutralino, which we assume represents the sole constituent of CDM in the 
Universe, by invoking a conventional Gaussian pdf in the likelihood (with an additional 
10% theoretical error). 



The constraints imposed by the recent LHC data are described in detail in Sec. 2.5 



2.5 LHC constraints 

The LHC at CERN started operations in September 2008, marking the beginning of an 
intensive period of investigation. The current plan (as of Summer 2011) is that the LHC will 
be running at the current centre of mass energy ^/s = 7TeV throughout 2012. Following 
this period, an approximately 18 month maintenance period is expected to commence in 
order to upgrade the LHC, after which the flrst attempts at collisions with ^/s = 14TeV 
are expected. 

Whilst the discovery of SUSY would obviously be of paramount importance, the task 
of actually identifying weakly-interacting massive particles (WIMPs) as the dominant con- 
stituent of DM using LHC data alone is challenging [^] , unless complementary information 



is provided by other experiments, e.g., by direct |[T^] or indirect [72] searches, or by speciflc 



accelerator signatures, such as the shape of the dilepton invariant mass spectrum in the 
MSSM focus point region |7|. 

To compute the region of the cMSSM parameter space where the LHC wih be able 
to achieve a 5a discovery, we consider studies conducted by the ATLAS collaboration for 
a centre of mass energy of 14TeV and an integrated luminosity (IL) of 1 fb~^ JT^ (which 
is expected to be achieved by 2014-2015) ^. We have adopted the analysis consisting in 
looking for events containing 4 jets + leptons -\-^t which have been performed for fixed 
Aq = 0,tan/3 = 10 values. Nevertheless these sort of analyses are relatively insensitive to 
tan/3 and Aq due to the fact that large tan/3 values, basically, alter the phenomenology 



of multilepton production channels [79] which are not considered here. We have followed 
the procedure described in [|^] for implementing the ^/s =14TeV and IL = 100 fb~^ LHC 
configuration (representing what might be achieved over a timescale of about 10 years). In 
this case, the authors have optimized the search looking for events containing a number of 
jets > 2 + n = 0, ..,6 leptons +^t where Aq = 0, tan/3 = 45 values have been assumed. 
Studies of this sort have shown that the detectability level exhibits only a mild dependence 



on Aq and tan/3 |81[]. Therefore, although the resulting LHC sensitivity is expressed as a 
detection region in the tuq, mi/2 plane for fixed values of Aq and tan/3, we assume it to be 
universal for all values of Aq and tan/3, as argued above. 



3. Results 

3.1 The cMSSM after null searches at the LHC 

We begin by displaying in Fig. |l| the favoured region of the cMSSM parameter space (in 
terms of the Bayesian posterior pdf) in the (ttiq, '^1/2) plane, where we have marginalized 
over all other parameters. The top and bottom rows show results for the flat and log priors 
respectively. From left to right, we have imposed current LHC constraints (left panels), 
and future constraints assuming that the LHC does not detect SUSY with ^/s =14TeV 
and IL^lfb"-*^ (middle panels) and ^/s =14TeV and IL^lOOfb"-*^ (right panels). This 
has been achieved by removing samples which lie below the 5a detection region, depicted 
by the magenta dot-dashed curve for the LHC with ^/s =14TeV and IL = lfb~^, and by 
the solid blue curve for ^/s =14TeV and IL= lOOfb"^ 

With the current LHC configuration, both the stau co-annihilation and /i-pole regions 
remain inaccessible. The latter (visible as a small, vertical sliver of probability probability 
at small values of mi/2) we see will become accessible with an integrated luminosity of 
lfb~^ and a centre of mass energy of 14TeV. For an LHC configuration of ^/s =14TeV 
and IL = 100 fb~^, we see that the stau co-annihilation region displayed becomes completely 
accessible. In case of a lack of detection, the surviving posterior probability is completely 



^We have also investigated the LHC configuration with ^/s =7GeV and IL = 1 fb ^, which is very close 
to what has already been achieved p&, iTA [r8|, and consequently, we omit our results for this case. 
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Figure 1: Posterior probability distribution for the cMSSM in the (mo, mi/2) plane, after null 
searches by the LHC with combinations of ^/s and integrated luminosities (IL) (left to right): 
Current LHC, 14TeV and lfb~^ and 14TeV and 100 fb~^, for the priors (top) and log priors 
(bottom). The encircled black cross represents the best-fit point. The inner and outer solid, black 
contours delimit the 68% C.L. and 95% C.L. posterior regions respectively. We also illustrate the 5cr 
detection threshold of the LHC by the magenta dot-dashed curve for ^/s =14TeV and IL = 1 fb~^, 
and by the solid blue curve for ^/s =14TeV and IL = 100 fb~^. 

confined to what, from here on, we call the funnel region, which encompasses the A- 
funnel and a small fraction of the focus point region which evades the current XENON-100 
constraints, since it is pushed to large tti^, for both choices of priors, as can be seen in the 
right-most panels of Fig. |l[ The percentage of the currently favoured cMSSM parameter 
space which would survive a lack of detection at the LHC is given in Table ^, for both choices 
of priors. Those values can be interpreted as probabilities for the nightmare scenario in the 
cMSSM: for the LHC configuration with ^/s = 14TeV and IL = 100fb"\ the probability 
of a non-detection lies between about 3% and 36%, depending on the choice of priors, with 
the flat prior giving a larger probability to the nightmare scenario (as a larger fraction of 
its posterior probability lies in the funnel region, beyond the reach of the LHC). 



3.2 Implications for Direct detection 

In this section, we investigate the impact of hypothetical null searches at the LHC for 
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Table 3: Percentage of the currently favoured parameter space in the cMSSM that would survive 
a lack of SUSY detection at the LHC with the given configuration of centre of mass energy ^/s and 
integrated luminosity IL. 
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Figure 2: Favoured region in the cMSSM once current constraints and future null searches at the 
LHC are taken into account, for flat priors (top) and log priors (bottom). We assume null searches 
at the LHC with the following combinations of ^/s and integrated luminosities (from left to right): 
Current LHC, 14 TeV and lfb~^, and 14 TeV and 100 fb~^. The encircled black cross represents 
the best-fit point. The inner and outer solid, black contours delimit the 68%C.L. and 95%C.L. 
posterior regions respectively. We also show the current 90% C.L. exclusion limit from XENON-100 
(magenta dotted), and the expected reach for for Phase 1 (solid red curve, expected to be reached 
by ^2012), Phase 2 (dashed green curve) and Phase 3 (dash-dotted cyan curve, expected to to be 
reached around 2020) future direct detection experiments. 



future ton-scale direct dark matter detectors. 

In order to establish the prospects for directly detecting cMSSM DM, we need to 
make some assumptions regarding the local neutralino density and velocity distribution. 
A consistent way to combine accelerator and DD data is to impose the 'scaling Ansatz' 



discussed in [|T5[ , which rescales the local density of neutralinos according to their predicted 
cosmological relic density. Here, however, we focus on the regions of the cMSSM parameter 
space where the neutralino is the sole constituent of DM, by imposing a Gaussian likelihood 
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Table 4: Fraction of the cMSSM parameter space surviving null searches at the LHC which can 
be probed by direct detection experiments or by IceCube (the latter after 5 years of observation). 
Those values represent the probability of discovery by the corresponding experiment (for each prior 
choice) assuming the LHC fails to discover SUSY at the given energy and integrated luminosity. 



on ft^h^ (see Table I), and therefore the scaling Ansatz becomes almost irrelevant. 

For an easier comparison with the existing literature we have fixed the astrophysical pa- 
rameters describing the density and velocity distribution of DM particles to the commonly 
adopted benchmark values: local CDM density P0,CDM = 0.4GeVcm~^; circular velocity 
vq = 235 km s"-*^ and escape velocity Vesc = 550 km s"-*^ (see, e.g., |2^ and references 
therein for a recent discussion of the astrophysical uncertainties relating to these quanti- 
ties). We have also neglected hadronic uncertainties in WIMP-nucleon couplings, adopting 
for the light quarks contribution to the nucleon form factors the values /tu — 0.02698, 
frd = 0.03906 and /ts = 0.36 [^. Accounting for such uncertainties by including them as 
nuisance parameters in the scan would only mildly change the numerical results presented 
here, but the main conclusions would remain unchanged, in agreement with the findings of 

In Fig.|2|, we display the results of the posterior probability distribution of the cMSSM 
in the (rn^, cr^p) plane following our scan corresponding to null searches at the LHC with 
the same priors and combinations of ^/s and integrated luminosities as displayed in Fig.[l|. 
For comparison, the expected sensitivity of upcoming direct detection experiments is also 
shown for the following experimental setups: 

• Phase 1: experiments will probe cross sections down to a^p ^ 10~^pb, correspond- 



ing to the projected sensitivity of SuperCDMS at SNOLab [|83| with a detector target 
mass of 27kg operating for 1 year (solid red curve), roughly equivalent to the reach 
of XENON-100 by the end of 2012. 



Phase 2: 



• if^nase z: experiments will probe cross sections down to cr^p ^ 10~^^ pb, correspond- 
ing to the projected sensitivity of SuperCDMS at SNOLab with a detector target 
mass of 145kg running for 3 years (dashed green curve). 

Phase 3: experiments will probe cross sections down to cr^p ^ 10~^^ pb, correspond- 
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ing to the projected sensitivity of SuperCDMS at SNOLab with a detector target mass 
of 1500 kg running for 4 years (dashed green curve), that should become available by 
2021 [|3[. The XenonlT is expected to reach a sensitivity of a^^ = 4 x 10"^^ pb 
by 2015 [84], and therefore to go down to cr^p ^ 10~^^pb on a similar timescale as 



SuperCDMS. 

The left panels of Fig. ^ show the favoured cMSSM region with current LHC data, 
for flat (top) and log (bottom) priors. Table|^ summarizes the percentage of our points 
that reside outside of the reach of the LHC, i.e., in the nightmare scenario (for a given 
configuration) that are detectable for each of the experimental direct detection phases 
described above. 

The /i-pole, stau co-annihilation and funnel regions displayed in the corresponding 
plots of Fig. H can be observed in Fig.|2| as "islands" in the parameter space bound by 
68%C.L. contours and spanning the approximate mass ranges: 50GeV < m^ < 60GeV, 
100 GeV <m^< 500 GeV and 200 GeV < m^ < 1 TeV respectively, with the best-fit point 
from current data (left panels) being located in the stau co-annihilation region for both 
priors. 

Large portions of these three regions are within reach of Phase 1 direct detection 
experiments, which will be able to detect between 30% and 68% (depending on the choice 
of priors) of the parameter space currently outside the reach of the LHC. For both choices 
of priors, we observe that Phase 2 experiments will cover a substantially larger fraction 
of the currently surviving parameter space, between 73% and 95% for flat and log priors, 
respectively. Phase 3 experiments are expected to cover entirely the most probable region 
of the cMSSM currently inaccessible to the LHC, a result that, as we shall see, holds true 
also for the ultimate reach of the LHC. 

After its ^18 month shutdown period starting in 2013, the LHC will be brought up 
to v^ = 14 TeV. The central panel of Fig. ^ shows the implications for direct detection 
assuming IL = lfb~^. The LHC data in this case will rule out the h-pole region and start 
to cut in to the stau co-annihilation region, pushing the neutralino mass to larger values 
(see Fig. 1^ below) and substantially reducing the percentage of surviving points within the 
reach of Phase 1 DD experiments, from 30% to 23% for flat priors, and from 68% to 35% for 
log priors. The prospects for detection with Phase 2 DD experiments remain fairly stable, 
while, again. Phase 3 DD experiments can probe all of the surviving parameter space. 



Lastly, for ^/s =14 TeV and IL = 100 fb~^, as expected from Fig. ^, for both flat and log 
prior scans the h-pole and stau co-annihilation regions are now fully accessible to the LHC, 
with the inaccessible part of the funnel region being quite similar for either set of priors. 
In both cases, we can see from Table^ that, whilst a small percentage (of order 1 — 2%) of 
points are accessible to the Phase 1 experiments. It is clear that Phase 2 detectors will be 
necessary in order to have a significant chance (i.e., of order ^ 30%, independently of the 
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Figure 3: Posterior probability distribution of the cMSSM, displayed in the (m^, Njf^) plane, 
corresponding to null searches by the LHC with combinations of ^/s and integrated luminosities 
(from left to right): Current LHC, 7TeV and lfb-\ 14TeV and lfb-\ and 14TeV and 100fb-\ 
when using flat priors (top) and log priors (bottom). The number of muon events N^S^ assumes 
5 years of observations by IceCube plus DeepCore. The encircled black cross represents the best- 
fit point. The inner and outer solid, black contours delimit the 68%C.L. and 95%C.L. posterior 
regions respectively. The horizontal dashed line gives the 90% C.L. detection threshold, and the 
encircled black cross is the best-fit in each scenario. 

choice of priors) of discovering dark matter in the context of the cMSSM, whose gaugino 
and scalar mass parameters would in this case be located in the high mass funnel region, 
if it evades discovery by the LHC in the long term. 

Interestingly, we find that even in the case of null searches at the LHC with ^/s =14 TeV 
and IL = 100fb-\ i.e., after many years of data taking after 2014, 100% of the cMSSM 
posterior will be probed by Phase 3 DD experiments, that will become available on a similar 
timescale. 



3.3 Implications for indirect detection with the IceCube neutrino telescope 

In this section we examine the prospects for detecting DM in the cMSSM with indirect 
searches in the nightmare scenario, focusing on the detection of high energy neutrinos from 
the Sun. Dark Matter particles are expected to scatter off nuclei in the Sun, and to sink at 
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Figure 4: ID pdf for the neutralino mass m^ for flat priors (top) and log priors (bottom), assuming 
Current LHC constraints (left), null searches at the LHC with y^ =14TeV and IL = 1 fb~^ (middle) 
and null searches at the LHC with ^/s =14TeV and IL = 100fb~^ (right panels). The encircled 
black cross represents the best-fit point. 

its centre, where they can subsequently annihilate. Among the annihilation products only 
neutrinos can escape, and produce a neutrino flux on Earth that is currently searched for 
with the IceCube neutrino telescope, located at the South Pole (see, e.g., [^]). IceCube 
is now fully built, with 80 strings in total, and it has been recently supplemented by a 
more densely instrumented region (consisting of six additional strings) at its core, called 
DeepCore (see, e.g., [|^, |8|]), which has the effect of lowering its neutrino energy threshold, 
E^^ to approximately 10 GeV and increasing the effective area for low-energy events [^]. 

The number of neutrino-induced muon events, Njp^ observable by IceCube from the 
direction of the Sun can be computed as 






(3.1) 



where ^^ is the incident muon neutrino flux at the detector, E^ is the energy threshold, 
and A^{Ejj) is the muon neutrino effective area of the detector for muon neutrinos with 
energy E^^ averaged over the northern hemisphere. The flux is then integrated over the solid 
angle AJ^q subtended by the Sun, corresponding to a cone with half-angle of approximately 
0.3°. 
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Figure 5: Summary of the reach of various probes in the m^ vs a^p plane for flat priors (left) 
and log priors (right) for samples resulting from a null search using the current configuration of 
the LHC. The solid, black contour delimits the 95% posterior region from current data (including 
Current LHC constraints). Equal weight samples from the posterior are shown and are coloured 
(and thinned by a factor of two for display purposes) as follows: green (magenta) squares are 
accessible to the LHC with ^/s = 14TeV and IL = lfb~^ (IL = lOOfb"^); blue squares are outside 
the reach of the LHC with ^/s — 14TeV and IL = 100 fb~^ (nightmare scenario); cyan filled circles 
are expected to be favoured at 90% C.L. by Planck; red crosses are accessible to IceCube. The 
yellow shaded region represents the reach of future Phase 3 direct detection experiments. 



The capture rate is computed using DarkSUSY v5 . 05 assuming, as above, a local CDM 
density p^ — 0.4GeVcm~^ and a Maxwell-Boltzmann velocity distribution with mean 
V — 235kms~^. In order to assess the prospects for discovering DM with IceCube, the 
expected signal, given by Eq. (|3.lD , must be compared with the the expected number of 
background events arising from the following sources: 



i) down-going muon events produced in a single cosmic ray shower in the southern hemi- 
sphere; 

ii) muon events, misreconstructed as upgoing, resulting from two simultaneous cosmic 
ray showers in the southern hemisphere; 

iii) an isotropic background of neutrino-induced muon events generated by atmospheric 
neutrinos; 

iv) neutrino-induced muon events generated by cosmic neutrinos; 
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v) and muon events generated by neutrinos produced by cosmic-ray interactions in the 
solar corona. 



We follow a methodology similar to that adopted in [jSg] and we only consider events 
from the direction of the Sun when it resides in the northern hemisphere, thereby eliminat- 
ing contribution from (i). The contribution from (ii) can be reduced to negligible levels by 
using appropriate analysis cuts and therefore it is reasonable to assume that IceCube can 
reject such a background [^]. The isotropic background contribution from cosmic neu- 
trinos (i.e., background (iv)) can be safely neglected (we estimate it to be approximately 
5 X 10~^ events per year), as well as the background contribution (v) from cosmic-ray in- 
teractions in the solar corona, which we calculate to be approximately 5 x 10~^ events per 
year. Our estimate of background due to up-going atmospheric neutrinos (iii) is obtained 



by substituting into Eq. (|3.lD the best-fit measurements for the isotropic atmospheric neu- 



trino flux from SuperKamiokande and AMANDA II presented in Fig. 10 of [Q, leading 
to an expected background, (n^), of approximately 4.9 events per year. From this back- 
ground, the number of muon events corresponding to a 90%C.L. detection threshold can 
be estimated using the Feldman-Cousins construction [^], giving Nqq = 32 events for a 5 
year period. 

In Fig. ^ we display the posterior probability distribution of the cMSSM, in the (m^, 
NJP) plane, corresponding to null searches by the LHC with combinations of y/s and 
integrated luminosities as in Figs.|| and |2|for flat (upper) and log (lower) priors on cMSSM 
input parameters. We also display the corresponding number of events associated with the 
90%C.L. Feldman-Cousins sensitivity estimate Nqq (dashed line). In TableQ we list the 
percentage of points inaccessible to the LHC that are detectable by IceCube (i.e., above the 
A^90 threshold) for each of our investigated LHC configurations. The left panels of Fig. ^ 
correspond to the current LHC configuration. For our log prior we can clearly identify the 
/i-pole, stau co-annihilation and funnel regions as the 68% C.L. regions spanning the mass 
ranges 50GeV < m^ < 60GeV, 100 GeV < m^ < 500 GeV and 200 GeV < m^ < ITeV 
respectively. The distribution of probability in these regions is extremely similar to that 
displayed in Fig. ^ for the (m^, a^p) plane. Despite the fact that each of these three regions 
are partially accessible to IceCube, only up to ^ 17% of the samples currently outside the 
reach of the LHC can actually be probed. For flat priors, both the stau co-annihilation 
and funnel regions are partially accessible (the /i-pole region is barely visible due to the 
resolution of the scan), but the relative suppression of the stau co-annihilation region with 
respect to the log priors leads to an even smaller percentage, ^ 10%, of detectable samples. 



The central panels of Fig.^ correspond to ^/s = 14TeV and IL = lfb~^. In analogy 
with the corresponding results in Fig.g, the /i-pole region and a large portion of the stau 
co-annihilation region are probed by the LHC, and therefore ruled out in case of null 
searches. This leads to a shift of the posterior to higher neutralino masses, and to an 
enhancement in correspondence of the funnel region, which tends to reside for linear priors 
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Figure 6: ID pdf for the lightest Higgs mass rrih for flat priors (top) and log priors (bottom), 
assuming Current LHC constraints (left), null searches at the LHC with ^/s =14 TeV and IL = 1 fb~^ 
(middle) and null searches at the LHC with ^/s =14 TeV and IL = 100fb~^ (right panels). The 
encircled black cross represents the best-fit point. 

at m^ ^ 500 GeV and Njp < 10, hence inaccessible to IceCube. For log priors, the 
neutralino mass is shifted towards lower values, and Njp towards slightly higher values, 
leading to a larger probability of detection, 11%, to be compared with 6.9% in the case of 
flat priors. 

Finally, for ^/s = 14 TeV and IL = 100fb~^ (right panels), IceCube will have little 
impact in probing the cMSSM in case of null searches at the LHC, for either sets of priors, 
resulting in detection probabilities in the range 0.1 — 1%. 

We have also verified that for all LHC configurations, the entire parameter space 
fraction accessible to IceCube will also be accessible to Phase 1 direct detection experiments 
(but not vice- versa). This means that IceCube on its own is not expected to be able to 
improve constraints with respect to even Phase 1 direct detectors. However, it will provide 
a valuable independent cross check in at least a subset of the parameter space. This 
is obviously important in terms of evaluating and reducing systematic effects, e.g., from 
astrophysical uncertainties, that would affect both direct and indirect detection limits. 

Several features in Figs. § and || can be better understood upon inspection of the ID 
pdf for 771^, shown in Fig. § for the same LHC configurations and same priors. The /i-pole, 
stau co-annihilation and funnel regions can be roughly identified in the left panels, and 
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the progressive shift towards large m^ as the LHC cuts into the stau co-annihilation and 
funnel regions is apparent from the central and right panels. It is also clear from this figure 
that the prior dependency of the results decreases as more and more constraining data are 
accumulated (left to right), as expected. 

Our results are summarized in Fig.||, which shows in a more concise way the param- 
eter space accessible to various observational channels in the (m^, a^p) plane, for both 
choices of priors. For this figure, we have additionally investigated which fraction of the 
parameter space that can be potentially ruled out by a more accurate determination of the 
relic abundance by the Planck satellite, which is expected to tighten current cosmological 
limits on the DM relic density by a factor of ^ 10 (see, e.g., [Q). We find that about 
95% of the cMSSM parameter space outside the reach of the LHC will be potentially ruled 
out by Planck (at the 90% C.L.), independently of the assumed LHC configuration. This 
is because models with the "correct" relic density (i.e., matching the current determina- 
tion by WMAP within Planck sensitivity) are essentially uniformly spread throughout the 
(7710,^^1/2) plane (see the cyan samples in Fig.|5|). Therefore, Planck has almost uniform 
power in probing them, independently of the parameter space fraction excluded by the 
LHC. Fig. P also clarifies that the reach of IceCube lies just below the current direct de- 
tection exclusion limits from XENON-100 (red crosses), while the LHC nightmare scenario 
(i.e., the region inaccessible to the LHC with ^/s = 14TeV and IL = 100fb~^, indicated 
here by the blue squares) is confined to higher neutralino mass regions of the parameter 
space, corresponding to spin-independent cross sections in the range a^p ^ 10~^ — 10~^^ 
pb. It is again apparent that Phase 3 DD experiments will be able to probe the entire 
surviving cMSSM parameter space, including the nightmare region. As mentioned above, 
further inclusion of astrophysical and hadronic uncertainties in the analysis is unlikely to 
change this conclusion, as the reach of Phase 3 DD experiments lies comfortably below the 
favoured parameter space. 

3.4 The role of the Higgs 

What is the role of the Higgs boson in this scenario? Obviously, a detection of the Higgs 
in the next few years would provide crucial information, but since we do not know if, and 
at which mass, it will be detected, we can only discuss what consequences the possible 
outcomes of Higgs searches would have on our scenario. 

To this aim, we have calculated the ID pdf for the lightest Higgs mass rrih, shown 
in Fig.^ for the same LHC configurations and same priors as in Figs.§ and ^. In the 
so-called decoupling limit, the lightest Higgs couplings to Standard Model particles are 
identical to those of the Standard Model Higgs [Q. This happens when uia ^ raz, a 
condition that is fulfilled for large values of mi/2- Upon inspection of Fig. ||, we see that 
the posterior is indeed pushed towards large values of 7711/2, especially in the right-most 
plots, corresponding to IL = 100 fb~^ of data at y/s =14TeV. 

A recent analysis of the prospects for detecting the Standard Model Higgs boson at 



ATLAS [94|, has shown that IL = 2 fb ■'^ of data at ^/s =8 TeV will be sufficient to exclude 
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Standard Model Higgs masses between 114 and 500 GeV at the 95% CL. Since our ID pdf s 
for null searches at the LHC with 14 TeV fall in the range 115-125 GeV, we can conclude 
that if the LHC will be able to rule out the Higgs in this range, then the cMSSM will be ex- 
cluded. On the other hand, a 5a discovery of the Higgs might take a much larger luminosity 
and might require running at larger energies. So if by the end of 2012 the LHC will have 
failed to rule out the last remaining window for the Higgs (currently, between 115 and 140 
GeV, considering the recently presented results from the ATLAS and CMS combination of 
1.6 + 1.6 fb~^ of data [^), then the future generation of direct detection experiment will 
have a good chance of discovering dark matter (within the cMSSM framework considered 
here). 



4. Discussion and Conclusions 



In this study we have quantitatively assessed the potential for discovering the cMSSM in 
case of null searches at the LHC (the "nightmare scenario" of particle physics) , with future 
direct detection experiments and with the IceCube neutrino telescope. 

Our main conclusion is that Phase 3 direct detection experiments (that can reach cross 
sections down to a^p ^ 10~^^ pb) will be able to probe entirely the favoured region of the 
cMSSM parameter space in the nightmare scenario of particle physics, therefore providing 
a unique opportunity to test SUSY even in case of null searches at the LHC. Interestingly, 
these experiments are expected to be built on a timescale of 5-10 years, similar to that 
required for the LHC to reach IL = 100 fb~^ of data at y/s =14 TeV. 

In our analysis, we fixed all astrophysical parameters to the standard values commonly 
adopted in the literature. In principle, however, one should take into account all particle 
physics and astrophysical uncertainties in order to combine in a self-consistent fashion 



accelerator and particle astrophysics experiments |T^, ^ Q. Including the uncertainty 
on the local DM density |9^, ^] would not change the fraction of cMSSM points accessible 
to direct detection experiments, since the SUSY points and the sensitivity curves would 
be rescaled by the same quantity (i.e., the ratio of the 'true' local density over the value 
adopted here). A more careful treatment of the velocity distribution and of the hadronic 
uncertainties in the neutralino-nucleon cross-section may instead have a stronger impact 
on direct and indirect searches [^, |9^, Q, but given the ample margin between the bulk 
of the nightmare cMSSM parameter space and the sensitivity of Phase 3 direct detection 
experiments, it is unlikely that our main conclusion would change significantly. 

Finally we have studied the implications of Higgs searches at the LHC, and we have 
shown that if a Standard Model-like Higgs is not found at the LHC with IL = 100 fb~^ of 
data at ^/s =14 TeV, that would basically rule out the cMSSM, while an actual detection 
in the appropriate mass range would provide additional motivation to continue the study 
of Supersymmetry with astroparticle experiments. 
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